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In aquatic plants such as the Characean algae, 
the force generation that drives cyclosis is lo- 
calized within the cytoplasm, yet produces fluid 
flows throughout the vacuole. For this to oc- 
cur the tonoplast must transmit hydrodynamic 
shear efficiently. Here, using magnetic resonance 
velocimetry, we present the flrst whole-cell mea- 
surements of the cross-sectional longitudinal ve- 
locity fleld in Chara corallina and show that it is 
in quantitative agreement with a recent theoret- 
ical analysis of rotational cytoplasmic streaming 
driven by bidirectional helical forcing in the cy- 
toplasm, with direct shear transmission by the 
tonoplast. 
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The pioneering work of Kamiya and Kuroda (Kamiya 
and Kuroda 1956) established that cytoplasmic stream- 
ing in Characean algae takes place not only in the cyto- 
plasm, but also produces fluid flows that extend through- 
out the entire vacuole. As pointed out long ago (Pickard 
1972) this result implies that the vacuolar membrane 
(tonoplast) efficiently transmits shear generated in the 
cytoplasm (Houtman et al. 2007) to the vacuole (Fig- 
ure la). While the ability of membranes to flow under 
shear is well-known from the "tank-treading" of erythro- 
cytes (Fischer et al. 1978), it is clear from recent stud- 
ies of lipid vesicles under extensional flows (Kantsler et 
al. 2007) that sheared membranes can undergo instabil- 
ities that lead to complex dynamics. With the ability 
to image the dynamics of the tonoplast directly in Ara- 
bidopsis by means of a green-fluorescent-protein-labelled 
tonoplast integral channel protein (Cutler et al. 2000), 
there is fairly clear evidence that such complex dynam- 
ics do take place in vivo^ and these may have implica- 
tions for flows throughout the cell. In Chara corallina^ 
for example, the actin- myosin system driving streaming 
is localized in two opposed helical bands at the cell pe- 
riphery, and the regions of high shear where they meet 
may create complex tonoplast dynamics. 



The possibility that streaming may impact on cel- 
lular metabolism by mixing the contents efficiently 
(Hochachka 1999, Pickard 2006) has been revisited 
recently through solutions of the coupled advection- 
diffusion dynamics of rotational cytoplasmic streaming 
(Goldstein et al. 2008, van de Meent et al. 2008). These 
calculations of velocities throughout the vacuole assumed 
that shear transmission by the tonoplast was complete, 
and showed that vacuolar flows can produce enhanced 
mixing, providing a possible mechanism for the home- 
ostatic role hypothesized by Hochachka. The first di- 
rect measurements of the wall-to-wall profile are those of 
Kamiya and Kuroda (1956), who studied rhizoid cells, 
'leaf cells sprouting from nodes, and internodal cells, 
and found a constant velocity within the cytoplasm and 
a curved shear profile within the vacuole. Mustacich 
and Ware (1976) improved on these measurements by 
using laser-Doppler scattering through a chloroplast-free 
window obtained by exposure to an argon laser prior to 
observation. Pickard (1972) obtained velocity measure- 
ments for a collection of native particles in an internodal 
cell of Chara braunii^ showing consistency with the ve- 
locity deduced under the approximation of non-helical 
indifferent zones. 




FIG. 1: Geometry of Chara corallina. (a) Cross-section of a 
Chara internode. Flow is driven in opposing directions along 
two helical bands separated by indifferent zones visible as light 
lines on the surface. Velocities are generated in a 10 /xm 
thick layer of cytoplasm at the periphery, but a shear flow 
extends into the central vacuole that takes up 95% of the 
volume of the cell, (b) Microscopy image of the sample prior 
to insertion into the glass tube. The dotted line shows a 
spiral of dimensionless wavelength X/R = 42. (c) Schematic 
of sample holder enclosed in a horizontal RF coil. The length 
of the tube is 40 mm. 
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FIG. 2: MRV measurements of the cross-sectional velocity profile of Chara corallina. (a) NMR image of the sample, with the 
measurement volumes at x = —3 mm , —1 mm, 1 mm marked in green, blue and magenta. Dotted lines mark the extrapolated 
position of the indifferent zones, as determined from the orientation of the velocity profiles at 5 measurement points between 
X = — 5 mm and x = 3 mm (marked with black arrowheads), (b) Cross-sectional plots of the velocity. The resolutions along 
the y— and axis are 16 /xm and 31 /xm respectively. Measurements were averaged along the x— axis over a Gaussian- weighted 
region of width 1 mm (full- width at half maximum). The fitted maximum velocities of the profiles are 45.2 /xm/s, 46.0 /xm/s 
and 46.9 /im/s. 



Here we report a technical advance in the study of 
cytoplasmic streaming by obtaining the fluid velocity 
throughout a cross-section of internodal cells in Chara 
corallina directly with magnetic resonance velocimetry 
(MRV). Our results allow for quantitative tests of re- 
cent fluid dynamical theories (Goldstein et al. 2008, van 
de Meent et al. 2008) and suggest further uses for mag- 
netic resonance imaging (MRI) in the study of large-scale 
streaming flows. 

In recent years, MRI techniques have increasingly 
found use in non-medical applications (Elkins and Al- 
ley 2007), and the advent of phase-shift MRV has made 
it possible to perform non-invasive flow measurements on 
microscopic scales (Callaghan 1994). In MRV the initial 
application of a pulsed magnetic field gradient encodes 
each spin with a 'label' describing its position along the 
direction of the gradient. At a time A later (the "obser- 
vation time"), a reversed gradient is applied, introducing 
a net phase shift in the orientation of the nuclear spin 
system that is directly related to the distance travelled 
in the direction of the gradient. By careful selection of 
the magnitude of the applied field gradients and the ob- 
servation time, velocities from ~ 10~^ — 10 m/s can be 
measured, depending on fluid properties. With sufficient 
time-averaging, spatial resolutions of 10 — 100 /im can 
be achieved, allowing imaging of biological systems on 
scales just slightly larger than those of typical single cells 
(Choma et al. 2006). While this technique has success- 



fully been applied at tissue level in a variety of plant 
systems (Scheenen et al. 2001, Kckenberger 2004, Windt 
et al. 2006), it is in the uniquely sized internodal cells of 
Chara that we can obtain measurements of flows internal 
to single cells. 

The Characean internode (Fig. la) is a single cylindri- 
cal cell with a diameter up to 1 mm and a length that 
can exceed 10 cm. The bulk of the volume of the cell is 
occupied by a large central vacuole, within which a 0.15 
M concentration of salts results in the 5 bars of osmotic 
pressure that lends the cell its rigidity. A layer of cyto- 
plasm ~ 10 /im in thickness encloses the vacuole at the 
cell periphery. Charophytes are recognised as the high- 
est genetic predecessors to vascular plants (Karol et al. 
2001) and in Chara most organelles common to higher 
plants are found in the cytoplasm. Somewhat uniquely, 
the millions of chloroplasts cover the cell wall, packed into 
helical rows that spiral along the inner surface (Fig. la- 
b). On the inside of those rows, bundled actin filaments 
act as tracks for myosins that drag structures within the 
cell (Kachar 1985, Kachar and Reese 1988) and thereby 
entrain cytoplasm. With streaming rates as high as 100 
/im/s, the myosin XI found in Chara is the fastest known 
(Shimmen and Yokota 2004). As a result of a reversed 
polarity of the actin filaments, the flow is organised in 
two opposing bands, producing a "barber-pole" velocity 
at the cell periphery. These bands are separated by indif- 
ferent zones, identifiable by the absence of chloroplasts 



3 



(a) (b) 




-400 -200 200 400 -400 -200 200 400 

x(|jm) x{|jm) 




-400 -200 200 400 -40 -20 20 40 

position (|jm) u^mr (Mm/s) 



FIG. 3: Comparison between measurements and theoretical solution, (a) Reconstructed profile created by aligning and 
averaging the data sets shown in Figure 2b. The black circle shows the position of the cell wall, (b) Corresponding theoretical 
profile, (c) Three wall-to-wall velocity sections at orientations of 0°, ±35° and ±70° with respect to the z— axis (see arrowheads 
in panels above). For non-zero angles the data is averaged over the section with positive and negative orientation. Markers 
show the measured values, lines show the theoretical prediction, (d) Comparison between measurements and theoretical values. 
Data points from all three profiles are shown as markers, with the color of the marker denoting the radial position of the data. 
For clarity, only 1 in 10 points is shown. 



and visible as two light lines crossing the surface of the 
cell. 

To measure the cross-sectional flow inside the intern- 
odal vacuole, a sample is inserted into a horizontal 
solenoidal radio-frequency (RF) coil (Fig Ic). Flow pro- 
files with a cross-sectional resolution of 16 /im x31 /im 
are obtained over a total period of 192 — 256 minutes, 
averaging over a volume of 1 mm along the longitudinal 
axis of the cell (see Materials and methods). Figure 2a 
shows the location of the five averaging domains on a 
particular sample, and Fig. 2b shows the velocity pro- 
files at three of those positions. Mean flow velocities are 
~ 45 /im/s, typical of a cell of radius R = 460 /im at 
the temperature of 298 K inside the coil (Pickard 1974). 
The helical pitch obtained from five short measurements 
2 mm apart (Figure 2a, arrowheads), is 18.6° /mm, giv- 
ing a dimensionless wavelength of XjR = 42.0, in good 
agreement with our estimate from the light-microscopy 
image in Figure lb. 

Aligning and averaging the three datasets in Figure 2b 
we further enhance the signal to noise ratio of the pro- 
file. The resulting cross-section (Figure 3a) is in excellent 
agreement with the hydrodynamic solution presented in 



earlier work (Goldstein et al 2008), shown in Figure 3b. 
In order to account for the ~ 20° of helical rotation along 
the averaging volume, we average the theoretical expres- 
sion (see Materials and methods) over a length of 3 mm 
with the same Gaussian weighting as the MRV measure- 
ments. The wall-to-wall 1-dimensional profiles in Figure 
3c show close agreement with the theoretical solution. In 
Figure 3d, measurement points are plotted against their 
theoretical values, in the manner of earlier velocity mea- 
surements (Pickard 1972). The color of the points indi- 
cates the radial distance from the centre of the cell. We 
see a remarkably good correspondence throughout the 
bulk of the cell, with deviations restricted primarily to 
points within a pixel from the cell wall, where partial- 
volume effects become significant. 

In summary, we have shown that magnetic resonance 
velocimetry can yield detailed information on the veloc- 
ity distribution of cytoplasmic streaming within single 
plant cells, allowing quantitative comparison with fluid 
dynamical theories. Natural extensions to this work in- 
clude studies of other streaming geometries found in na- 
ture and of the spread of tracers as a probe of mixing 
(Esseling-Ozdoba et al. 2008). 
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Materials and methods 

Plant materials 

Chara corallina v. australis was obtained from the 
Botanic Garden of the University of Cambridge, cour- 
tesy of J. Banfield, and grown in a non-axenic culture, 
rooted in non-fertihzed soil in a 100 liter tank filled with 
Artificial Pond Water (1 mM NaCl, 0.4 mM KCl, 0.1 
mM CaCl2). The tank was kept at room temperature 
and illuminated with a bench lamp on a 16/8 hour day- 
night cycle. During illumination, the light intensity at 
the top of the tank was ~ 250 lux. Samples of suitable 
size were placed in a Petri dish under a microscope to 
verify healthy streaming. 

In preparation for measurements, the sample was in- 
serted into a 3.0 mm outer diameter, 1.6 mm inner diam- 
eter glass capillary 40 mm in length, which was pre- filled 
with Forsberg medium (Forsberg 1965). The capillary 
tube was closed with PDMS plugs and a small volume of 
silicone grease to ensure a good seal. 

Magnetic resonance velocimetry 

Velocimetry experiments were performed on a Bruker 
Spectrospin DMX 200, 4.7 T magnet with a 20 mm long 
solenoid coil of diameter 3 mm. images were acquired 
at 199.7 MHz. Spatial resolution was achieved using 3- 
axis shielded gradient coils providing a maximum gra- 
dient strength of 49 G cm~-^ in each direction. Trans- 
port is measured over the observation time, A, and since 
the RMS displacement increases as A^-^ due to diffusion 
and as A^ due to convection, for short A diffusive (inco- 
herent) displacements can dominate over convective (co- 
herent) ones, particularly in slowly convecting systems. 
Therefore, to weight the measurement towards the con- 
vective field a stimulated echo sequence was used to en- 
able a large observation time for motion encoding. Fur- 
ther spatial imaging gradients were applied after the mo- 
tion encoding to minimize diffusive attenuation. Hard 
90° RF pulses were used except for the final pulse which 
was a Gaussian-shaped selective 90° RF pulse 512 /is in 
duration. Experimental parameters used for the veloc- 
ity images were: observation time, A = 500 ms; veloc- 
ity gradient duration, d = 1.62 ms; gradient increment, 
^inc = 10 G cm~^; number of velocity gradient incre- 
ments, 2; recycle time, TR = 1.9 s; number of scans 
= 16; field-of-view = 2 mm x 2 mm; pixel array size 
A^read X ^phase = 128 X 64; iu-plauc Spatial resolution = 
16 fim X 31 /im; slice thickness = 1 mm; measurement 
duration = 64 min. 

For each data set in Figure 2b, 4 measurements were 
taken between 3 and 7 hours apart and subsequently av- 
eraged. A Gaussian smoothing with 31 /im FWHM was 
selectively applied to the streaming region whilst avoid- 
ing blurring at the cell periphery. The streaming ve- 
locity in the individual measurements was not observed 
to vary significantly over the 15 hour period of acqui- 
sition. Throughout the measurements the temperature 



was maintained at 298 K by a regulated airflow system. 

Hydrodynamic solution 

The hydrodynamic solution was obtained as a mode 
expansion of the equations for Stokes flow (Goldstein et 
al. 2008). The system of equations assumes an infinitely 
long helical cell, allowing simplification of the problem to 
a 3-dimensional flow with symmetry, where the velocity 
depends only on the radial coordinate r and a helical 
angle (p = 6 — {27r/X)x. Here, is the polar angle in the 
7/z-plane and A is the wavelength of the helical bands. 

The solution of the flow problem can be expressed as 
a sum over a series of modes /f(r)sin((2z + 1)^), with 
the radial modes fi{r) given in terms of a sum of Bessel 
functions. In the results presented here, 64 modes were 
used for the expansion. The velocity at the boundary was 
taken piece- wise constant on the bands, with a smooth 
crossover at the indifferent zones given by tanh((/?/e) and 
tanh((7r — (p)/e) respectively. The cross-over width used 
was e = 11 jam. 

To allow better comparison with the MRV measure- 
ments, a series of profiles was averaged along the x-axis 
over a length of 3 mm, using a Gaussian weighting with 
a spread of 1 mm full-width at half maximum. 
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